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Abstract 

Spontaneous emission of an atom (molecule) placed near a nanocylin- 
der of elliptical cross-section of an arbitrary composition is studied. The 
analytical expressions have been obtained for the radiative and nonradia- 
tive channels of spontaneous decay and investigated in details. 



1 Introduction 

As was first pointed out by Purcell pQ one can change the atomic decay rate by 
using a resonant cavity. It is well known now that not only resonant cavities, 
but also any material body may exert influence on the spontaneous decay rate 

Moreover, control for the spontaneous decay rate is widely used in practice 
at elaboration of new efficient light sources [3] . 

Due to rapid development of nanotcchnologics influence of nanoparticles 
of different shapes on radiation of atoms, molecules and nanocrystal quantum 
dots is of key importance. This direction often referred to as nanophotonics 
0]. Influence of spherical [5] [7], spheroidal [5] -[111] and even ellipsoidal [IT] 
nanoparticle on radiation of atoms and molecules is investigated in details. 

Of special interest is the influence of dielectric fibers or metallic wires of 
different cross-sections on the decay rate of a single atom. This is because the 
charged wires were employed successfully in order to control the atomic motion 
[P2l [13] . The cylindrical geometry is also important for the investigation of 
fluorescence of substances in submicron capillaries [141 115j . A very important 
application of the decay rate theory in the presence of dielectric fiber had been 
in the case of the photonic wire lasers [Jj)]-[TH]- Finally, the circular cylindrical 
geometry appears naturally at analysis of the carbon nanotubes [50] [UJ . 

The influence of a perfectly conducting cylindrical surface on decay rates 
has been well investigated for an atom placed both inside a cylindrical cavity 
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[22]- [22] and near a perfectly conducting cylinder [53]. Spontaneous emission of 
an atom located inside a coaxial waveguide was considered in |27j . 

The interaction of an atom with a dielectric, semiconductor, or a metallic 
cylinder is a more complicated process. Van der Waals interaction between an 
atom and a circular metallic nanowire was considered in [25]. The first inves- 
tigation of the decay rate of an atom placed at the axis of a circular dielectric 
fiber was undertaken within a classical approach in [29] [30]. Recently that prob- 
lem attracted new interest [3T]-[53]. In [371 [35] the spontaneous emission near 
carbon nanotubes was considered. 

Recently non-circular nanocylinders (the so called nanobelts or nanoribbons) 
were successfully synthesized [3S]-[13]. Figure 1 illustrates the experimentally 
obtained nanobelts [39] . 

The significance of studying such nanoobjects is very great because they 
form the basis for the synthesis of all kinds of nanodevices, such as nano-scale 
transducers, actuators, or sensors [44]- [47]. The influence of such structures on 
optical characteristics of atoms and molecules had not been yet considered, as 
far as we know. 

The aim of this work is to investigate theoretically the influence of dielectric 
and metallic cylinders of elliptical cross-section on the spontaneous emission of 
atoms and molecules placed close to the cylinder. This geometry is close to 
synthesized nanobelts and nanoribbons. 

Though it is possible to separate variables in the Hclmholtz equation within 
the elliptic cylinder coordinates, still the electrodynamics and optics of the di- 
electric elliptic cylinder remain insufficiently studied, due to mathematical dif- 
ficulties concerned with the Mathieu function [48 . In this paper we shall try 
to make up for this deficiency, and shall report the investigation results deal- 
ing with the spontaneous emission of the atom placed near an infinite elliptic 
cylinder. 

Some elements of the theory for a spontaneous decay of an atom near ma- 
terial bodies will be presented in section 2, where the main attention will be 
paid to the nanobodies, whose dimensions are small as compared to the radia- 
tion wavelength. In section 3, we report solution on a problem of spontaneous 
decay of an atom near the perfectly conducting elliptic cylinder with arbitrarily 
sized cross-section when retardation effects can be important. In section 4, the 
radiative and nonradiative channels of a decaying atom located near an elliptic 
nanocylinder, which is made of an arbitrary material, will be considered analyt- 
ically within the quasistatic approach. The results obtained in sections 3 and 4 
will be analyzed and illustrated graphically in section 5. The geometry of the 
problem under study is illustrated in Fig. 2. 
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2 Some elements of the theory for the sponta- 
neous emission of an atom placed near the 
nanobodies 

At weak coupling of an atom and a nanobody, i.e. in the case of the exponential 
spontaneous decay, the expression for the lincwidth 7 a of excited level a has the 
form [35] : 

2 3 

7a = 7a,o + ^ E E d MV m {GMr\ r'; u na )} 6 (co an ) (1) 

n a.fi=l 

where 7 0j o is the linewidth in free space without nanobodies; the matrix 

element of the dipole moment operator between the states a and n; oj an — 
(W a — W n ) /Ti] and G aj 3, the reflected part of the classical Green's function 
connected with the reflected field E r a of a dipole dg a is placed at point r' by the 
relation 

3 

^(r)=E^G a/3 (r,r> na ) (2) 

(3=1 

In Eqs. fl} and ([2|), a, (3 — 1, 2, 3 are the indices of the Cartesian axes, and O is 
the Heaviside's unit-step function. Note that Eq.fl]) has been obtained within 
the framework of the most general assumptions (Fermi's golden rule and the 
fluctuation-dissipation theorem), and has a very large domain of applicability. 
From Eq.JT]) it is seen that total linewidth is the sum of partial widths, and 
below we shall consider, for simplicity, the rate of transitions between some two 
states only. 

In the case of the nanobodies one can use the Rayleigh's perturbation long 
wavelength theory. The Green's function of the reflected field may be repre- 
sented as the power series in the wave vector |50] : 

G aP (r, t';cj) = (r , r') + kG% (r, r') + fc 2 G$ (r, r') + tk 3 G^ (r, r') + . . . (3) 

where G^l (j—0, 1,2,...) are the coefficients that can be found from a solution 
of the corresponding quasistatic problems; k — lo/c. The first three terms ([3]) 
describe the near fields, while the radiation fields emerge starting from the 
4 th term. By substituting Eq.Q into (JlJ for the total rate of spontaneous 

decay 7 e tg from the state e into g near the nanobody, we obtain the following 

expression: 
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7e^g = T 




} 



nonradiativc 



} (4) 



a,f3=l 



radiative 



Here A: » wo/c, and 70 and o^o are unperturbed by a nanobody the rate and 
frequency of transitions between the states a and g. The first term of Eq.(TJ} is 
nonzero for the absorbing media only, whereas the rest of the terms are nonzero 
in the absence of absorption as well (the dielectric or the perfectly conducting 
nanobodies). Those terms describe mostly the radiative losses. Thus, to deter- 
mine the main terms of the nonradiative and radiative losses it is necessary to 
determine G ®1 (r',r') and G^l (r',r'), respectively. 

It is a complex problem to determine directly the radiative losses described 
by the 3-d order terms by k. But in the case of an atom localized close to 
a nanobody the radiation is dipole-like, and the total dipole moment of the 
atom + nanobody system could be found from G^l (r,r') at large distances 
from the system. Thus, in the case of the nanobodies the radiative rate of the 
spontaneous transition (the transition from e into g) will be described by [35] : 



where d to tai = do + Sd is the total dipole moment of the atom + nanobody 
system, and Sd is the dipole moment induced in the nanobody. Here we will 
omit superscript eg in transition dipole moment d^ 9 and in 7 e yg and 70, e >9 . 

It should be noted that beside the radiative and nonradiative chanels of 
the spontaneous decay for a dielectric cylinder (and for any dielectric waveg- 
uided structure) one must take into account energy transfer into the waveg- 
uided modes. In case of the dielectric (nonmetallic) nanobodies that correction 
is exponentially small [35] and may be neglected. But in case of the complex 
dielectric permittivity of a nanocylinder (real metals), the spontaneous decay 
waveguiding rate may be no longer a small value, and that correction should be 
taken into account. 

Thus, to describe the decay rate of an atom in the presence of a nanobody 
one should determine the reflected field E' r in a quasistatic approximation, and 
the dipole moment <5d induced in the nanobody. For a nanocylinder which size 
is comparable with radiation wavelength one should use full set of Maxwell's 
equations to find the reflected field. 



radiative 



|d + M| 2 

4 



(5) 
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3 Spontaneous decay rate of an atom near a per- 
fectly conducting elliptic cylinder 

The problem of the spontaneous atomic decay near a perfectly conducting ellip- 
tic cylinder of any geometry may be solved analytically in the elliptic cylinder 
coordinates (Fig. 3) (1 < u < oo, — 1 < v < 1, — oo < z < oo) [3FJ 

x = fuv, y = fy (u 2 — 1) (1 — v 2 ), z — z, (6) 

where / = a 2 — b 2 is the half-distance between the cylinder foci. Expression 
u = uq = a I f specifies the surface of the studied cylinder. 

Total field from the arbitrary current j (r) and charge p (r) densities near the 
perfectly conducting elliptic cylinder can be found from the equations [511 152) 

Ej M (r) = J dr'[ik j -^l+p{T')-^G™{v,v'), 

H T Z E (r) = \ j dr' [j(r') x V] z G TE (r,r') , (7) 

where the electric and magnetic Green's functions G™ and G TE are of the 
following form [51] (u > u' > Uq ) 



er=e,on=0 " ^ ■> > 

xH<U* n (pf,u) (Ja n (0f,u') - H^a n 09/, u') 

J(T n ((3f,U ) \ 



G TE (r/) = 



V V [ dae ia ( z - z ') S(Jn ^ S(Tn 
^ ^ J Ma n (8f) 

xiJ«<7„ 09/, u) (ja n 09/, t*') - ff (1 V„ 09/, «') 



x 



(d/dw)HWa n (0f,w) 



(9) 



1 /2 

where /3 = (fc 2 — a 2 ) ; Sct„ denote the "angular" even (cr = e) and odd 
[a = o) Mathieu functions defined as follows [48] : 
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Se 2n (h, cos 0) = J2 Be 2n ( h ) cos ( 2m9 ) ' 

m=Q 

oo 

Se 2n +i (h, cos 6) = Be 2n+i ( h ) cos (( 2 ™ + !) ) . 

771=0 
OO 

So 2n+1 (ft,cos0) = J] Bo^-H W sin ((2m+l) 0) , 

m— 

oo 

So 2 „ (ft, COS 0) = B °2n CO Sin ( 2m61 ) ' ( 10 ) 

where Ba™ (ft) are the expansion coefficients (Boq 711 (ft) = 0); Ja n and H^a n 
are the "radial" Mathieu functions, which are expressed through the same co- 
efficients Bcr™ (ft). For example, functions Ja n can be written down as [38] : 

— oo 

^£(-lf~ m 5e^ (ft)J 2m (ft M ), 

— oo 

m=0 
/7T v'?-! 



Jo 2 „ + i (ft,u) 



2 u 

£ (-1)-™ (2rn + 1) B<%$t (ft) J 2m+1 (hu) , 



x 

m=0 

— / 2 _ 1 °° 

£ (-l)" _m 2mBo^ (ft) J 2m (fttt) , (11) 

m=l 

where J n are the Bessel functions. For the functions H^'a n the analogous 
expansions will be valid provided that the Hankel functions are used (substitute 
J m — ► H$ in Eq. pip ). The normalization constant Ma n (ft) is determined by 

the expression Ma n (ft) = J dd [Sa n (ft, cos 6)] . 

o 

Integration over the longitudinal wave number a in Eqs. ([8J and ([9j) is 
performed along the real axis assuming that the wave number k has a positive 
infinitely small imaginary addition (k — > fc + i0 + ). It should also be noted 
that the first part of Eqs. (J5J and (which is proportional to Ja n (//?, u')) 
corresponds to Green's function for free space (in the absence of a cylinder), 
whereas the second part (which is proportional to H^a n (//3, u')) corresponds 
to the reflected field. 

Note also that for a near-circular cylinder (/ — > 0) one can obtain from Eqs. 
(J5J) and ([9]) the respective expressions for the Green's functions in cylindrical 
coordinates (p, ip) if we set that 



G 



Se n {f(3,v) i fn s cos (rup) 



VMe n (fp) V 2 (1 + 8o,n) 

So n (f(3,v) sin (nip) 

--Jo n (ff3,u) -> -=— Jn(pp), (12) 



\fMo~jJ&~) W V2 

where i5o,n is Kronecker's <5-symbol. Analogous expressions for the functions 
H^Un are obtained from (jT3J) via a substitution J n — > if^. 

The charge and current densities at excitation by a point dipole (atom) have 
the form: 



p = - (d V) £(r-r') e- lwt , 

j = -^ a(r-r')e- lui , (13) 

where do is the transition dipole moment (below just the dipole moment) of an 
atom placed at point r' (below, the primed coordinates will denote the atomic 
position). The exponential term e~ luJt will be omitted. 

Knowing the field components ([7]) we shall determine the rest of the field 
components both for TM and TE modes [S3J 



E™ = —I— , / : t^±- E 1 M = 

^u,a Q1f\ „.2 „.2 ' '-'v.a 02/ 



lOi / U 



2 -l dE™ TM _ i a I 1 - v 2 dE z a 



P 2 f V u 2 — v 2 du v ' a (3 2 f V u 2 — v 2 dv 



„, , jfc / 1 - v 2 dH^ E j-,TE ik I u 2 -1 dH^ E 



E 1E = —\- - E 1 * = ——\— - — (14) 

I3 2 f\lu 2 -v 2 8v ' v ' a f3 2 f\u 2 -v 2 du ' y ' 

In Eq. (|14p the a index denotes the respective Fourier transform for the z coor- 
dinate. Total field is the sum of the TE and TM field components. 

In the case of the dipole source (fT5)l . the expressions J7]) for the fields may 
be rewritten in the following form: 



E™ (r) = fc 2 d , z G™ (r, r') + (do V) ^,G™ (r, r') , 

Hj E (r) = -ik[d Q x V'] 2 G T£ (r,r'). (15) 

By using (HJ) , ([9]) , and (fT4f , (fT5)) one can represent the spontaneous emission 
rate (Q} as the sum of the TM and TE contributions: 

^ = 1 - ^Re (W™ + W TE ) . (16) 
7o 2 

In some particular cases Eq. fTo]) will take the next relatively simple form. 
Dipole moment is oriented along x-axis (do = doe x ) 
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A) Atom is placed on the x-axis (r' = x'e x ) 



W 



TM 



oo J. 

4^ / da 

n=o{ 



k 3 (pjf Me n (f3f) 



_d 

du 



7 H^e n ((3f,u') 



-i 2 



Je n ((3f, u ) 
HWe n ((3f,u )' 



(17) 



W. 



TE 



da 



u ' 2 - 1 ^{ k(f3f) 2 Mo n ((3f) 



(d/dw) Jo n (/3f,w) 



(d/dw)HWo n (J3f,w) 
B) Atom is placed on the y-axis (r' = y'e y ) 



(18) 



W. 



TM 



4 

772 



V / da- 



a 2 l(d/dv>)Se 2n+1 (Pf,v>)\ v , =0 ] 



Me 2n+1 (0f) 

2 Je 2n+1 (J3f, u ) 
# (1) e 2 „+i (/3/,«o) 



oo K 
n=0 

x [i/ (1) o 2 „ (/3/X)] 2 



a 2 [(5/di/)So 2 „(/3/y)U c 



Mo 2 „ (/?/) 



■/Q2n (/?/, ^p) 

HWo 2n ((3f,u Q )' 



(19) 



da [Se 2n (pf,0)Y 
k(Pff Me 2n (Pf) 



— H We 2n W,u') 



-i 2 



(d/dw) Je 2n (0f,w) 



+4 1 



(d/dw) HWe 2n (f3f,w) 
da [So 2n+1 (/?/,0)] 2 



W—Uq 



d 



g^HW 02n+1 (f3f,u>) 



hi kWf M °2n + lW) 

2 (d/dw) Jo 2n+1 (f3f,w) 



(d/dw)HWo 2n+1 (/3f,w) 
Dipole moment is oriented along y-axis (dp = dpe y ) 



(20) 
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A) Atom is placed on the x-axis (r' = x'e x ) 



W 



TM 



oo " 

da 

n=0 J n 



H^o n (J3f,u') 



£3 tfff Mo n [fi '/) 
J On (/?/, Up) 



(21) 



TE 



. v - ^ f da 1 

h{ kWfMe n W) 



(d/dw) Je„ (f3f,w 



(d/dw)H^e n (0f,w) 
B) Atom is placed on the y-axis (r' = y'e y ) 



(22) 



W, 



TM 



4 1 



-' 2 )S/ 



« 2 [Se 2w (/3/,0)] : 
fc3(/?/) 2 Me 2li (/?/) 



Je2n (/?/, Uo) 



#«e 2 „ (/3/,n ) 

+4 ( 1 _±)y[ da gg [g^gM: 



du 



7 H^o 2n+1 (pf,u') 



J02n+1 (fif. Up) 

HWo 2n+1 (f3f,u ) : 



(23) 



TV, 



TE 



,'2 



oo " 



d« [{d/dv')Se 2n+1 {(3f,v')\ v , =0 } 



„=o S kWY 



H^e 2n+1 {(3f,u') 



Me 2n+1 {(3 f) 
2 {d/dw) Je 2n+1 {(3f,w) 



{d/dw)H^e 2n+1 {(3f,w) 



-V 



da [{d/dv')So 2n (Pf,v')\ v >= 



„_„■< Hdff Mo 2n {f3f) 



H^o 2n (J3f,uj 



{d/dw) Jo 2n {(3f,w) 



{d/dw)HWo 2n {f3f,w) 



(24) 
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Dipole moment is oriented along z-axis (do = d$e z ) 
A) Atom is placed on the x-axis (r' = x'e x ) 



W. 



TM 



n=0 



P 2 



k 3 Me n 08/) 
B) Atom is placed on the y-axis (r' = y'e y 



Je n {(if, u ) 
H^e n (pf,u o y 



wT E = 0. 



(25) 



w: 



TM 



W. 



TE 



n=0 



da 



(3 2 [Se 2n ((3f,0)} 2 
k 3 Me 2n (/3f) 



H {1) e 2n (j3f,u') 



Jd2n (/?/, U ) 



+4^ da 



n=0 



P 2 [So 2n+1 (/3/,0)] : 



fc 3 Mo 



2n+l 



H^e 2n (fi f, uo) 
H^o 2n+1 (pf, u')Y 



Jo 2n +i (flf. Up) 
H^o 2n+1 (/3f,UoY 



0. 



(26) 



Case of low ellipticity ( f — > 



If a cylinder has low ellipticity (/ — > 0) one can find explicit expressions 
for the spontaneous emission rate (|16p in the form of the power series /. That 
can be done if the expansion for the coefficients Be™ (h) at h — ► is known. 
Keeping the terms of the second-order smallness by h on can write down the 
expressions for -Bcr™ (/i— > 0)in the following form (n — 1, 2, 3, . . .): 



Be 



Be" 



2(ra-l) 



Be 



2(n+l) 



Be 2 



/i 2 



16 (2n - 1) ' 
fr 2 

16(2n + l) : 



Be 



2n 



1 - 



h 2 



8(4n 2 -l)' 



(27) 



R e 2»-1 
-° e 2n+l 



Be 



2n+3 
2n+l 



ft 2 , 

1+ 32' 

/7 2 

— Be 2n+l 
32n' g2 " +1 



32 (n+ 1)' 



32' 

1 - 



32n (n + 1)' 



(28) 
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and 



Bo\ 



r>_2ra— 1 
B °2n+1 



3/t 2 
"32"' 

h 2 



Bo{ ss - 



32 : 



32n (2n+l) 



So: 



2n+l 



1 



2n+1 ~ 2n + 1 



h 2 



32n(n + 1) 



Bo: 



2n+3 
'2n+l 



/, 2 



Bo 



2(n-l) 

In 



32 (n + 1) (2n + 1) 
// 2 



(29) 



Bo, 



2(n+l) 
2n 



32n (2n- 1)' 
h 2 



Bo 2n 



1 

2n 



8 (4n 2 - 1) 



32n (2n+l) ' ^ 

The rest of the coefficients Be™ are either zero or have a higher infinitesimal 
order. By substituting Eqs. (f27|) - (|30|) into (fT0|) and (jTTJ) , and assuming that in 



the cylindrical coordinates (p, 99) we have /« w p ^1 + | sin 2 93 J and w w 

cos </? ^1 — 5 ^) sm2 i we can nn d expansions in the region of / — > for the 

Mathicu functions and their derivatives with accuracy up to f 2 . By restricting 
ourselves to the main expansion terms only, and assuming that fu' = p' and 
fuo = a we obtain the following expressions for W™ and W TE . 
A) Dipole moment is oriented along p-coordinate line 



W, 



TM 



00 K . 

(2-*o,n) / da 

n ** 



W, 



TE 



n=0 



a 



dz 



HP (*) 



z=f3p> 



J n (f3a) 
((3a) 



o U 2 ) • 



da 



(d/dz) J n (z) 



(d/dz) (z) 



k {(3 P >y 

B) Dipole moment is oriented along (^-coordinate line 

2 J n (Pa) 



z—(3a 



(31) 



TM 



2 > n~ 

71=1 



da- 



F (/3p') 

A' 



hP W) 



w, 



TE 



«=0 ^ * 



— (a 



i?r\ 1} ((3a) 



1 2 



+ o(/ 2 ), 



z=0p' 



(d/dz) J n (z) 



(d/dz) H n 1] (z) 



0(f 



(32) 



z—f3a 
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C) Dipole moment is oriented along z-axis 



k 

32 



n=o q (paj 

T^J B = 0. (33) 

Note that the main terms (|3~lj) - ([3"3"]) coincide as in the case of the analogous 
expressions for W p , W v , and W z in the circular cylinder of a radius 26J. This 
circumstance supports the validity of our calculations. 

Atom is placed in close proximity to surface of a cylinder 

In the most interesting case of an atom placed in close proximity to the 
surface of a cylinder (u' — > uo ). The decay rate of a tangentially to surface 
oriented dipole tends to zero. The nonzero rate of decays is possible if the dipole 
moment orientation is normal to the cylindrical surface, i.e. for the dipole that 
is oriented along the u-coordinate line 



W 



TM 



,'2 



EE 

a—e,o n— 



da- 



a 



[Sa n {Pf,v')Y 



d 



H^a n ((3f,u') 



TE 



= 4 



du 
1 - v 12 



-i 2 



U—Uq 

k 



EE 

<r=e,o n=0 



Ma n {fif) 

Ja n (Pf,u ) 
H^o-n 08/, uo)' 

da [{d/dv')Sa n W,v')} 2 



H^a n (/3f,uo) 



k{fif) 2 Ma n (f3f) 
2 (d/dw)Ja n (0f,w) 



(d/dw) H^a n (Pf,w) 



(34) 



When elliptic cylinder radius is tending to zero kfuo — > 0, the decay rate (|34p 
of a tangentially to surface oriented dipole tends to zero the main contribution 
to the decay rates of normally oriented dipole is due to the terms with n = 
and n = 1 of (|34[) . By expanding the expressions (TIT))) and pT|) into series at 
kf, kfuo — > (the procedure of finding the expansions for the coefficients -Be™ 
and £<7j™ is described in detail in [15]) one can obtain from (JTHJ) and (|34|) the 
next simple asymptotic formula: 



12 



1_ 

,7o, 



kfu' —kfuQ — >-0 



2(u 



(kfY 



1 H — arctanL* (uq) 

7T 



(u + \/ul - l) 



,12 



4 (In 2- 1) 
7r 2 (l + £* 2 M) 



(35) 



in which 



L* (uo) = 



In 



4 



«o-l 



"7-E 



(36) 



where -fE ~ 0.57722 is the Euler's constant. In a particular case of a circular 
cylinder of a-radius (/ — > 0, fuo = a) we obtain from (|35|) the expansion that 
coincides with analogous asymptotic expression for the circular cylinder [26j 



kp' —ka — >0 



2 (ka)" 



1 H — arctanL* 

7T 



4 (In 2- 1) 

~TT 2 (1+L* 2 ) 



(37) 



where L* = f [ln(f ) + 7£ ] . 

In the most interesting case of a very thin nanobelt (uq — > 1) the decay rate 
of the atom (having normal orientation of the dipole moment) can tend both 
to infinity (the atom is at the nanobelt edge (v' = ±1)) and take finite values 
at other positions of the atom. This signifies nontrivial character of the found 
expressions ([35]). 



4 Spontaneous decay rate of an atom near el- 
liptic nanocylinder made of arbitrary material 

In the case of the dielectric or metallic cylinder of arbitrary dimensions and 
arbitrary permittivity e the problem of the spontaneous decay rate of an atom 
may also be solved analytically [54]. But such a solution is so cumbersome 
that it is very difficult to use it practically. In this section we restrict ourselves 
to a practically important case of a nanofiber or a metal nanowire of elliptical 
cross-section. It may be easy to describe the spontaneous decay, in this case. 

4.1 The rate of the radiative decays 

In order to find the spontaneous radiative decay rate of an atom placed at point 
r' near to elliptic nanocylinder of arbitrary material it is necessary to calculate a 
dipole moment induced in a nanocylinder, in accordance with Eq.([5]). In order to 
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find the dipole moment one may use the found analytical solution [TT] for three- 
axis ellipsoid and turn one of the semi-axes to the infinity (without using the 
elliptic cylinder coordinates). As a result, the expression for the induced dipole 
moment on the elliptic cylinder made of material with the dielectric constant e 
with semi-axes a and b on the x- and y-axes, respectively (Fig. 2), will take the 
form: 



3 
2 

d , x x' d 0y y' 



Sd = -- (d , x P xx F a (£') + d ^P yy F b (£') e y 



vV + e) (6 2 + V) w + e & 2 + C 

P xx x' P xx y' \ ( x' 2 y' 2 Y 1 . . 



a 2 +e b 2 + i> v V( fl2 +^') 2 (& 2 + e) 5 

where the components of the polarizability tensor of an elliptic nanocylindcr 
are: 



P xx = \ab(a + b) 



3 a + sb' 

1 w .N e-1 

-a& (a + 6 — , 

3 ea + o 



P yy = ^abia + b)——^, (39) 



where £' is the positive root of equation a ? + g, + gjf^ = 1; 

MO : 2 



^(O = ; 7 (40) 

One can obtain the following expressions for some particular cases of the 
dipole moment orientation along the Cartesian axes. 

A) Dipole moment is oriented along x-axis (do = doe x ) 



. , do (e — 1) ab (a + b)e 

da = — 



V^^Tf (V a2 + I' + V&Tf) {a + eb) 
do (e - 1) ab {a + b) x' ( x' 2 y' 2 



(a 2 + ef 2 v/^^+e V (° 2 + t'f ( &2 + O s 

(a 2 + e) (a + eb) + (b 2 + ?) (ea + b) 



(41) 
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B) Dipole moment is oriented along y-axis (do = doe y ) 



<Sd = - 



d (e — 1) ab (a + b) e y 



V^Tf (Va 2 +¥+ V^Te) (ea + b) 
d (e — 1) ab (a + b)y' I x' : 



+ 



(a 2 + Y) (a + eb) + (b + ?) {ea + b) 
C) Dipole moment is oriented along z-axis (do — doe z ) 



(42) 



Sd = (43) 

Note that the expressions (|4"Tj) and (|42|) in case of circular nanocylindcr 
(a = b) are significantly simplified and take the well known form [35] 



Sd 



e- 1 
e+l 



^) [-d + 2 (d n') n'] 



(44) 



where d = d , x e x + d 0:V e y ; n' = (x'e x +y'e y ) / p' and p' = y/ x' 2 + y' 2 . 

Because the induced dipole moment (|4ip - ((^5|) is known, one can find the 
radiative rate of the spontaneous decay by using Eq.©. Let an atom be placed 
on a;-axis at point x' . Then it is sufficient to consider two orientations of the 
dipole moment, in accordance with Eqs. (j4Tj) - (|4"3)) . 

A) Dipole moment is oriented along x-axis 



radiative 



ab (a + b) 



£-1 



y/x' 2 - a 2 + b 2 (x' + \]x' 2 - a 2 + b 2 ) \a + eb 
B) Dipole moment is oriented along y-axis 



(45) 



radiative 



l - 



ab (a + b) 



\Jx' 2 - a 2 + b 2 (V + \Jx' 2 - a 2 + b 2 ) \ea + 



e - 1 



(46) 



If the dipole moment of an atom is oriented along the nanocylinder axis (z- 
axis) then the radiative rate of the spontaneous decay should not be changed. 

If an atom is inside a dielectric elliptic cylinder without losses, then the 
expression for the radiative rate of spontaneous decays will take the following 
form at arbitrary dipole moment orientation, do = do. x e x + do. y e y + do. z e z , 



radiative 



Eb 



2 ^2 



l 0, X 
~dj 



ea + b 



2 J2 



d 2 
"n 



d 2 

I 2- ' 



(47) 
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which is independent of the atomic position inside the nanocylinder. In deriving 
(f47|) we neglect the local field factor [55]. It should be noted that formula ([47)) 
in agreement with that obtained in [56j for the case of two-dimensional ellipse 
(to obtain that result we must substitute in (|4"T)) d . z = 0). It also be noted, 
that expression (ITT)) can be applied to describe of an atom in elliptic cylindrical 
cavity in material with permittivity e$. In this case one should make substitution 
e — > 1/eo and multiply the obtained expression by ^/eq. 

Note that Eqs. (|45)) and (|46| have been derived as a limiting case of a 
finite-size ellipsoid with fixed permittivity. These formulae are unsuitable for 
describing the radiative decay rate of atoms near conducting (metallic) cylinder, 
because the waveguided modes (plasmon modes) of the infinite cylinder had 
not been taken into account. The correct expression for the case of perfectly 
conducting cylinder is given in Eq. (|35p . 



4.2 The rate of nonradiative decays 

In accordance with section 2, the rate of nonradiative decays of an atom near a 
nanobody is determined by a quasistatic solution of the problem near an elliptic 
nanocylinder. A Green's function of the Laplace equation expressed in ellipti- 
cal coordinates may be obtained from the Green's function of the Hclmholtz 
equation (the first terms in the right-hand side of Eqs. ((8)) and ()9))) in the limit 
k — > 0. As a result of such a procedure in the region uq < u < u' the expression 
for the scalar Green's function may be written in the form 



Go (r,r') 



r - r'| 

OO 

Sa n (iaf, v) Sa n (iaf, v') 



X 



4z / da cos (a (z — z')) 

a—e,o n—0 



Ma n (iaf) 

xJa n (iaf,u)H ( - 1) a n (iaf,u') (48) 

It should be noted here that function Sa n (iaf, v) , as follows from (flO)) and 
(fTT)) . is the real function because the coefficients Ba n (iaf) are real. Func- 
tions J o n (iaf,u) and H^'a n (iaf,u) are either real or imaginary depending 
on the parity of n. The product Ja n (iaf, u) H^a n (iaf, u') is always purely 
imaginary. Hence, the right-hand side of Eq. (|4"8")) is the real function. 

Using Eq. (|4"8")) one can obtain a quasistatic solution for the potential of a di- 
electric nanocylinder, with dielectric constant e, excited by the point-like dipolc 
source located at the point with (u',v',z) coordinates. The solution for the 
potential <p r induced in space near the cylinder has the form: 



OO 

at a- fA y?>\ ST ST f j ii i\\ S(Jn (iuf,v)S(Jn(iaf,v') 
= MdoV) ^ daco S (a(z-z)) — 

(J — C j O Tl — q 

xH {1) a n (iaf, u) H {1) a n (iaf, u) Aa n (iaf, u ) , (49) 
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where the elliptical wave reflection coefficient has the form 



d 

Aa n (s,u ) = (e - 1) — Ja n (s,w) 

ow 



d_ 

dw 

d_ 

dw 



H^a n (s,w) 
Ja n (s,w) 



Ja n (s,u ) 



Ja n (s,uq) 



H^a n (s,u ) 



(50) 



W — Uq 



For metallic elliptic cylinders with low losses (e = e' + is"; e' < Q, e" <C |e'|) 
denominator of (|50|) becomes close to zero for some longitudinal wavenumbers 
(plasmon modes). In this case quasistatic approximation is not valid and one 
should use full set of Maxwell's equations. In the case of elliptic nanocylinder 
made of real metals (e = e' + is", e" ~ 1) or dielectrics the denominator of 



(I50p is always nonzero and quasistatic approximation is valid in this case. So 
we restrict ourselves for the case of dielectric elliptic nanocylinders made of real 
metals or dielectrics. 

By calculating the reflected field E r = — Vcjf with the help of and using 
([2]) and Qwe obtain the following expressions for the quasistatic contribution 
to the nonradiative rate of the spontaneous decay of an atom near the dielectric 
(metallic) nanocylinder. 

Dipole moment is oriented along x-axis (do = doe x ) 



nonrad 



k 3 P 



EE ^ 



a—e,o n— 



,'2 



[Sa n (iaf, v')} 



Al2 



Ma n (iaf) 

(d/dw)HM(r n (iaf,w)\ 
HWa n (iaf, u') 

1-v' 2 \ (d/dw)Sa n (iaf,w)\ w=v , 



H^o-n (iaf,u') 



So- n {iaf, v') 
xRe [Aa n (iaf, u )} . 

Dipole moment is oriented along y-axis (do — do e y ) 



(51) 



6( u '2_x) (1-v' 2 ) 

k 3f2 ( u ,2 _ „/2)2 



EE 

a— e,o n— 



da 



[Sa n (iaf,v')Y 
Ma n (iaf) 



H^o-n (iaf,vf) 

, (d/dw) Sa n (iaf, w)\ 
Sa n (iaf, v') 



2 f , (d/dw) H^a n (iaf, w)\ 



HWa n (iaf, u>) 

2 

Re[Aa n (iaf,u )}. (52) 
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Dipole moment is oriented along z-axis ( do = d$e z ) 



\ nonraa n uxj . r n ,. „ .-.w 

7 1 6 / daa 2 l Sa n [laf, V )\ 



7o / , k 3 t—' ^— ; 7 M(T„ (ior/ ) 

,u/ z <r=e,ori=0 n 



2 [5ct„ (to/, u')]" 



fe 3 ^ d «« M«r n (ia/) ■ (53) 



If an atom is placed very close to the nanocylinder the nonradiative spon- 
taneous decay rate may be approximated by the expressions, which describe 
nonradiative processes near a plane surface [58] . For dipole with tangential and 
normal to surface of cylinder orientations we have respectively 



) nonrad , . / -. 

tang 16(fcA') 3 U+l 

\ nonrad / \ nonrad 

7 \ o / 7 



^ 70 / norm \ 70 / tang 

Here A' is the distance from atom to surface; in the case of an elliptic cylinder 
we can write A' = C \J ^r+^ jx + T^r^p'; where C' is the positive root of the 
equation + = 1. 

5 Graphic illustrations and discussion of results 

Let us first consider an atom near a dielectric nanocylinder. Figure 4 illustrates 
the radiative rate of the spontaneous decay (|4"5|) and (|4^|) near the dielectric 
nanocylinder. As an atom moves away from the surface the decay rate tends 
to the value, which is corresponds to the case of an atom in free space in the 
absence of a nanobody. As the dielectric constant of the nanocylinder material is 
increased, the radiative rate for the normally oriented dipole moment of an atom 
increases by tending to large (finite) value, and for the tangentially oriented 
dipole moment it decreases by tending to low values (zero). The most increase 
in the decay rate occurs for the normally-oriented atom in the direction of the 
larger semi-axis (Fig. 4(a)). 

From Fig. 5 it is seen how the decay rate changes for an atom located in 
close proximity to the surface of a nanocylinder depending on the ratio between 
the a and b elliptic semi-axes. If one of the nanocylinder semi-axes tends to 
zero (the case of a dielectric nanobelt) then from (jl5|) and (jlp| it follows that 
^radiative _^ £ 2 f QJ . fog normally-oriented dipole moment of an atom placed 
at the belt edge, while for the tangential orientation of the dipole moment we 
obtain Y adiative /~fo -> 1/e 2 (Fig.5(a)). For the atom placed at the belt facet, 
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independent of the dipole moment orientation, we obtain ^radiative _^ ^ 
(Fig.5(b)). 

Figure 6 illustrates radiative and nonradiative rates of the spontaneous decay 
as the function of the distance of an atom to the surface of the dielectric elliptic 
nanocylinder having the permittivity e = 3 + ilO~ 8 . The atom is shifted along 
the x-axis and its dipole moment is oriented along this axis. To calculate the 
Mathieu functions we have used of the method [59] based on the solution of the 
generalized eigen- value problem for the coefficients Bo™. The integration of 
(j5"T]) - (f53"]) was performed by the mean- value method using the Aitken's process 
[BU] for acceleration of the integral convergence. As seen from Fig. 6 the radiative 
rate of the spontaneous decay for a dielectric with a very small imaginary part 
of permittivity is dominating, as a whole, over a nonradiative rate, and takes 
maximal values at the cylinder surface. But if the atom is placed very close to 
the nanocylinder surface, the nonradiative rate exceeds over the radiative one, 
and at v! — uq it turns to the infinity, as follows from the asymptotic expression 
(j5"4"]l depicted by dotted line in Fig. 6. Let us determine the distance, on which 
the radiative and nonradiative rates for the normally oriented dipole moment of 
an atom coincide. By equating (|45|) with (|54|) . for the case of a small imaginary 
part of the permittivity e = e' + is" of a cylinder (e" « e'), we get that at 




(55) 

the nonradiative contribution to decay rate is equal to the radiative one. For 
the case represented in Fig. 6 we obtain: (x' — a) /a as 4.4 x 1CP 3 , which means 
that one may actually neglect here the nonradiative processes. If the atom is 
y-shifted and has the dipole moment oriented along this axis, then from (|5"5"]) we 
obtain 




(e'a + bf 



4(e' + l) e' 2 (a + by 



1/3 

. . (56) 



Figure 7 illustrates decay rates for a silver nanocylinder with e = —3.02 + 
i0.21 (A = 373 nm [61]). The ratio of the elliptical semi-axes is a/b=3 and 
a=6 nm (ka as 0.1). As is well seen, the region nearby cylinder, where the 
nonradiative rate is greater than radiative, is essentially wider that the region 
represented in Fig. 6. This is due to a larger imaginary part of permittivity 
of the nanocylinder material. At small distances from the cylinder surface the 
nonradiative rate of the spontaneous decay increases in accordance with (|54| 
(dotted line of Fig. 7), and at large distances from the cylinder surface, it tends 
to zero. 

Now let us consider the case of a perfectly conducting cylinder. The integra- 
tion of Eqs. (fT7"]) - ([2"6"]) was performed by the method of mean values with using 
the Aitken's process. Figure 8 illustrates the spontaneous decay rate of an atom 
placed near a perfectly conducting elliptic cylinder (ka = 1), as the function 
of the atom's position along the larger semi-axis a. As is seen, the decay rate 



19 



achieves the maximum in close proximity to the cylinder surface (v! — uq) for 
the normally oriented dipole moment of an atom (x-line of Fig. 7(a) or y-line of 
Fig. 8(b)). For the tangentially oriented dipole moment we obtain zero at the 
surface of the cylinder (y-line of Fig. 8(a) or x-line of Fig. 8(b), and also z-line) 
resulting from a full compensation of the dipole moment of an atom. As the 
ratio of semi-axes increases (when the smaller semi-axis b decreases) the max- 
imum rate of the spontaneous decay is increased (cf. x- lines in Fig. 8 (a) and 
Fig. 8(b)) and tends to infinitely large values in the limit of b — ► 0, which agrees 
with the asymptotic formula (|35p for the nanocylinder. 

Figure 9 illustrates the spontaneous decay rate of an atom placed near a 
perfectly conducting elliptic cylinder (ka — 1) as the function of the atomic 
displacement along y-axis (towards a smaller semi-axis) . The decay rate achieves 
the maximal value (normal orientation of the dipole moment) and the minimal 
value (tangential orientation) at the surface of a cylinder. At high values of 
the semi-axes ratio, the curves y and z (Fig. 9(b) for a/b — 100 3> 1) look like 
the well known solutions [58j for the spontaneous decay rate of an atom near 
a perfectly conducting plane surface with the normal (curve y) and tangential 
(curve z) orientations of the dipole moment. The maximum value of the decay 
rate here is about 2. The behavior of decay rate of an atom with x orientation of 
dipole momentum is slightly different from the case of the plane mirror because 
of edge effects in case of strip. 

Figure 10 illustrates the decay rate of an atom placed in close proximity to 
the surface of a perfectly conducting cylinder (u' = uq). The dipole moment is 
oriented normally to the cylinder surface along the larger semi-axis (Fig. 10(a)) 
and along the smaller semi-axis (Fig. 10(b)). As seen from Fig. 10, the expression 
(|35j) obtained in the case of ka — > for the atom with normal orientation, tends 
to infinitely large values, which disagrees with the finite values derived from the 
quasistatic expression (|45[) in the limit of e — > oo. This difference follows from 
the fact that the amplitude of the surface current propagating along the cylinder 
at ka — + increases infinitely large [2"5] . In other words, perfectly conducting 
cylinder of infinitely small cross-section is a kind of antenna, which is effectively 
excited by a point dipole (atom). This means that the induced dipole moment of 
the perfectly conducting cylinder and the rate of spontaneous decay of an atom 
also tend to infinity at ka — ► (see Eq. ([36|) ). The analogous process is also 
observed in the case of a circular nanocylinder (dotted line in Fig. 10, Eq.((37|)). 

Figure 11 illustrates a comparison for the rate of spontaneous decay f(T6]) 
obtained by numerical integration (I34|) and an asymptotic expression (|35| for 
the atom with the normally oriented dipole moment in close proximity to the 
surface (u 1 = uq) of perfectly conducting cylinder. As it seen, the asymptotic 
expression agrees well with the numerical data up to ka — 1. It can, therefore, be 
suitable for a fast evaluation of the spontaneous decay rate near the nanocylinder 
and for verification of the numerical solution. 
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6 Conclusion 



Thus, in this paper we obtained the analytical expressions for the radiative 
rate of the spontaneous decay of an atom near a dielectric nanocylindcr of 
elliptical cross-section within the framework of a quasistatic approximation. 
The retardation effects are taken into account for the spontaneous decay of an 
atom near a perfectly conducting elliptic cylinder, where the analytical solution 
can be found at arbitrary parameters of the problem. Special attention is paid 
to the case of a nanocylinder where the largest semi-axis of the cylinder is 
much less than wavelength. It is shown that a decrease in the cross-section of a 
metallic (perfectly conducting) nanocylinder leads to an infinite increase in the 
spontaneous decay rate for a normally oriented dipole moment of the atom. A 
simple asymptotic expression for the spontaneous decay rate has been obtained 
in the case of a perfectly conducting nanocylinder. For the atom placed in close 
proximity to a belt edge one can obtain arbitrary large values of the decay rate 
by means of tuning elliptic semi-axes. The obtained results can be used for 
checking the numerical solution in a more complex case of the dielectric elliptic 
cylinder. The results of this work may be also helpful at studying an interaction 
of atoms (molecules) with nanobelts and nanoribbons. 
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List of Figure Captions 



Fig.l The nanobelts [39]. 

Fig. 2 Geometry of the problem. 

Fig. 3 Elliptic cylinder coordinate system. 

Fig. 4 Radiative decay rate of an atom near a dielectric elliptic nanocylinder 
as the function of the position of an atom, (a) atom is shifted along x-axis; (b) 
atom is shifted along y-axis. Dipole moment orientation: x - along x-axis; y - 
along y-axis. Solid line corresponds to e = 3; dotted line corresponds to e = 10. 
Nanocylinder semi- axes ratio a/b = 5. 

Fig. 5 Radiative decay rate of an atom near the surface of a dielectric (e = 3) 
elliptic nanocylinder as the function of the semi-axes ratio for the atom in close 
proximity to the surface of the nanocylinder. (a) atom is placed on x-axis at 
point x' — a; (b) atom is placed on y-axis at point y' — a. Dipole moment 
orientation: x - along x-axis; y - along y-axis. 

Fig. 6 Radiative and nonradiative decay rates of an atom as the functions 
of its position relative to a dielectric elliptic nanocylinder with permittivity 
e = 3 + il0~ 8 . The ratio of the nanocylinder semi-axes a/b — 6; ka = 0.1. The 
atom is shifted by x-axis with the dipole moment oriented along this axis. The 
dotted line corresponds to the asymptotic form (|54p. 

Fig. 7 Radiative and nonradiative decay rates of an atom as the functions 
of its position relative to the silver elliptic nanocylinder with permittivity e — 
—3.02 + i0.21 (A = 373 nm [61]). The atom is shifted by x-axis with the dipole 
moment oriented along this axis. The ratio of the nanocylinder semi-axes is 
a/b = 3; a = 6 nm. The dotted line corresponds to the asymptotic form ([54]) . 

Fig. 8 Decay rate of an atom placed near a perfectly conducting elliptic cylin- 
der (ka — 1) as the function of atom's position relative to the cylinder. The 
atom is shifted along the x-axis. Dipole moment orientation: x - along x-axis; y 

- along y-axis; z - along z-axis. The semi-axes ratio: (a) a/b = 3; (b) a/b = 10. 

Fig. 9 Decay rate of an atom placed near a perfectly conducting elliptic cylin- 
der (ka = 1) as the function of atom's position relative to the cylinder. The 
atom is shifted along the y-axis. Dipole moment orientation: x - along x-axis; y 

- along y-axis; z - along z-axis. The semi-axes ratio: (a) a/b = 3; (b) a/b = 100. 

Fig. 10 Decay rate of an atom placed in close proximity to the surface of a 
perfectly conducting elliptic cylinder (u' = uq ), as the function of ka. The ratio 
of the cylinder semi-axes: 1 - a/b = 3; 2 - a/b = 10; 3 - a/b — 100. The dotted 
line corresponds to the case of a/b — 1. (a) atom is located on x-axis at point 
x = x', and its dipole moment is oriented along x-axis; (b) atom is located on 
y-axis at point y = y', and its dipole moment is oriented along y-axis. 

Fig. 11 A comparison of an asymptotic expression for the spontaneous decay 
rate (|35|) with (fT6|) obtained by numerical integration of (f34l) for the atom with 
the normal dipole moment orientation placed in close proximity to the perfectly 
conducting cylinder surface (u' = Uo )• A solid line corresponds to a numerical 
solution; the crosses denote asymptotic expression. The ratio of semi-axes: 1 - 
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a/b = 3 and 2 - a/b = 1. Position of an atom relative to the cylinder surface is 
marked by arrow. 
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